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Abstract 

We consider exact solutions generated by the inverse scattering tech- 
nique, also known as the soliton transformation. In particular, we study 
the class of simple real pole solutions. For quite some time, those solutions 
have been considered interesting as models of cosmological shock waves. 
A coordinate singularity on the wave fronts was removed by a transforma- 
tion which induces a null fluid with negative energy density on the wave 
front. This null fluid is usually seen as another coordinate artifact, since 
there seems to be a general belief that that this kind of solution can be 
seen as the real pole limit of the smooth solution generated with a pair 
of complex conjugate poles in the transformation. We perform this limit 
explicitly, and find that the belief is unfounded: two coalescing complex 
conjugate poles cannot yield a solution with one real pole. Instead, the 
two complex conjugate poles go to a different limit, what we call a "pole 
on a pole". The limiting procedure is not unique; it is sensitive to how 
quickly some parameters approach zero. We also show that there exists 
no improved coordinate transformation which would remove the negative 
energy density. We conclude that negative energy is an intrinsic part of 
this class of solutions. 
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1 Introduction 



1.1 Background 

We consider exact solutions generated by the inverse scattering technique, also 
known as the soliton transformation, in the so-called cosmological case (see 
section ^| below) . We do not consider the stationary axisymmetric case since we 
are interested in solutions of wave-like character and cosmological application. 
A comprehensive review of this subject and applications related to that of this 
paper has been given by Verdaguer Q. In the review one also finds references 
to earlier work in the field by Ernst, Chandrasekhar and others in the 1970s 
and 1980s. A general discussion of waves modelled by exact solutions is given 
in the book by Griffiths 

The inverse scattering technique was first applied to general relativity in 
1978 by Belinskii and Zakharov ||. Their algorithm generates exact vacuum 
solutions to Einstein's equations from a known "seed solution" , also called the 
"background" . This background is multiplied by the so-called scattering matrix 
to yield the new metric, called the perturbation. Note that the use of the 
word perturbation does not imply an approximation, the new metric is still 
an exact solution. Any poles in the scattering matrix give rise to solitonic 
solutions, and these have earned the inverse scattering technique its other name 
"soliton transformation" . The transformation is closely related to the Backlund 
transformations, and the exact relation has been clarified by Cosgrove[|J |5). 

1.2 Simple real pole solutions 

The inverse scattering technique can recreate well-known solutions. One inter- 
esting example is the Kerr-NUT metric generated from flat space [[| with two 
real poles in the scattering matrix. Homogenous cosmologies of Bianchi types 
I through VII can be generated with the inverse scattering technique. New so- 
lutions have also been produced, often with gravitational solitons. One class of 
such solutions that has been studied extensively |l], [| 0, ||, |j| is the "one real 
pole" perturbation on a Bianchi I or II background (for the purposes of this 
article, we choose the Kasner vacuum metric as a background). One real pole 
solutions may have merits as cosmological models; to see this we briefly survey 
the physical behavior of the solutions. 

The relevant region for cosmology is t > (see fig. |l|). A solitary shock wave 
emerges from the initial singularity at z — wq, t = 0. Behind the wavefront, the 
shock wave leaves way for the background (Kasner) metric. 

Thus this soliton-like disturbance sweeps away some initial disorder which 
surrounds the singularity. This has earned the single real pole solution generated 
from a Bianchi I or II metric its nickname "the cosmic broom" . (From the Italian 
la scopa cosmica.) Simple real pole solutions generated from other seed metrics 
or with more poles display similar behavior. The cosmic broom represents an 
interesting exact-solution model of gravitational waves of cosmological origin. A 
background of such waves could in principle be detected, and the cosmic broom 
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Figure 1: The light cone structure. Pole at wq. 



would provide an intriguing classical alternative to waves generated from the 
initial conditions of the universe in quantum cosmology. 

Furthermore, it turns out that the perturbation is undefined in region I, so 
we match the perturbation in region II across the light cone to the background 
(Kasner) metric in region I. While other choices for the matching are possible, we 
match to the background for two main reasons. One reason is that the perturbed 
metric in region II itself goes into the unperturbed one when approaching the 
light cone. The other reason is that we wish to maintain the traveling wave 
interpretation, which is most meaningful when the metric is asymptotically equal 
to the background on both sides of the disturbance. 

In general, one real pole is not sufficient to recover the background at spatial 
infinity. Instead, this can be achieved with two or more distinct real poles Q. 
Since each real pole defines a light cone, one then has to perform a number 
of matchings. Some of these matchings will involve matching the background 
within the inner light-cone to one-pole solutions. Consequently, a first step is 
to match the background to a one-pole solution, even though one is ultimately 
interested in n-pole solutions or a stochastic background of waves. Therefore, 
the problem described in this paper can be seen as part of the more complicated 
matching procedure for several distinct real poles. 

A good motivation to use a Kasner-type background is that it has been 
shown that any "generic" singularity is followed by a succession of Kasner eras 
[ jTo| . Moreover, a Bianchi I metric is approximately unaffected close to the 
singularity by the presence of matter JTT) . (Recently, doubt has been cast on 
this somewhat folkloristic claim Jl2|.) 
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1.3 Problems 



The cosmic broom suffers from metric coefficient infinities on the light cone in 
canonical (Belinskii-Zakharov) coordinates. Recently it was shown that these 
coordinate infinities can be removed by a certain (singular) coordinate transfor- 
mation (|, IJ. However, a new but potentially less serious problem was created 
by this transformation. Although the metric is now continuous, the Ricci tensor 
goes distributional-valued, i.e. acquires 5-functions, at the light cone. This kind 
of (5-function can normally be interpreted as an impulsive wave, e.g. a null fluid, 
(see the book by Griffiths ^ on the subject). Unfortunately, the null fluid in 
this case would have negative energy density ||. Let us exhibit this explicitly 
for the left hand side of the light cone [z < wq), a similar discussion holds for 
the right hand side. From Einstein's equation with the aforementioned Ricci 
tensor one finds 

T V = ~7T~l — 5(t + z — wo) Vi vj 

OTTt g t t 

where v a — l/^/2g u (d/dt — d/dz) is along the light cone. An observer at 
constant z < wq with 4- velocity £ a = (1/ -Jgtt) d/dt would measure the following 
energy density when the wave front passes: 

P =T ij ee = -—±—5{t+z- WQ ) < o 

167rf g tt 

which is manifestly negative and thus violates the weak energy condition. (Our 

signature is (H ) as is usual for inverse scattering applications, so gtt > 0. 

Also we are at t > 0.) Therefore, one would naively think that the present 
"null fluid" is merely an artifact of the singular transformation and that the 
"actual" one real pole solution has no null fluid. (See also work by Gleiser and 
Dfaz|13 14], [l5j on the problem of matching solutions in different regions and for 



methods of removing the singularities on the light-cone and obtaining smooth 
extensions.) 

There are at least two ways to find a more physical form of the one real pole 
solutions. The most obvious is to "sew together" the coordinate patches in a 
different way, i.e. to find a better transformation which is non-singular but still 
removes the coordinate infinites which the metric has in canonical coordinates. 
The coordinate transformation is discussed in section ^. 

As an alternative, there seems to be a general belief ([[0, sec. 3.7.1) that 
the one real pole solution be sought as the real limit of two complex conjugate 
poles wq ± it. The complex conjugate pole solution is known to be smooth 
everywhere. Indeed, this limiting procedure has often been implicitly invoked 
to support the validity of simple real pole solutions such as the cosmic broom. 
(The reason one does not simply use the complex pole solution directly is that 
it contains disturbances which travel at v > c, so it in itself does not admit a 
traveling wave interpretation similar to that of the real pole solutions, where 
the shock wave is strictly located to the light cone and thus travels at v = c.) 
The limiting procedure is the subject of section ||. 

Our conclusions are that 
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• The limiting procedure actually yields a different solution, so it cannot be 
used to explain away the null fluid, and 

• There is no improved transformation, so this avenue out is also closed. 

2 Setup 

We use the conventions and notation of Belinskii and Francaviglia |7|. The 
inverse scattering technique is applicable for vacuum metrics of the form 

ds 2 = f{z,t)(dt 2 - dz 2 ) - g ab (z,t)dx a dx b 7 a, b =1,2, (1) 

which describes a spacetime admitting an Abelian G 2 group of isometries ||. 
Examples of such metrics are given in the introduction. This form of the metric 
is in the "cosmological case" . The other (stationary axisymmetric) case, which 
we do not consider, is related to this one by a complex transformation, but 
the stationary axisymmetric solutions are quite different in character. In the 
following g is the 2x2 matrix representation of g ab - 

Belinskii and Francaviglia |Q found an explicit expression for g and / of 
the complex conjugate pole solution. Note that we have replaced a = t in 
their expression. Also note that the superscripts "(2)" and "(0)" refer to "two 
poles" and "no poles" (background), respectively. A bar represents complex 
conjugation. 



= 5 i° fc ) + ^(M-M)(t 2 -|M| 2 ) 



-(t 2 - fi 2 )Q 2 2L a L b 



\(t 2 ~ [i 2 )QnL a L b - (fj, - n)Q22Qiig^! 
A* 



(2) 



m = C9 m (3) 

3 2 (t 2 -[ifi) 2 (t 2 -n 2 )(t 2 -fi 2 )(n-fi) 2 t 2j K ) 

where the vector L a and the scalars Qki are 

L-x — mx9n ' ^2 — m 2 g22 , 

n — m 2 n^ _1_ m 2 (11) 
Wll — m x9xX + m 2.9 2 2 ) 

r\ (°) 1 (0) r\ 

Q 12 = rrumxg^' + m 2 m 2 3 22 ; = Q21 , 

Q22 = fh\gf^ + m\g { 22 = Qxx > 

the mi and m 2 are determined by the seed metric (s is the Kasner parameter): 

m a = ( m i(2wp) _s , m Q2 (2wfi) s - 1 ) (4) 

with mob arbitrary complex constants, and T> is found from 

A = Q11Q22 - Qu — t 2 (miffi2 - wim 2 ) 2 

V = t 2 Q 11 Q 2 2(fi-fi) 2 + A(t 2 -^ 2 )(t 2 -fl 2 ) . (5) 
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The pole "trajectory" fi(z,t), called trajectory because it depends on z and t, 
is 

fj,(z, t)=W- Z± y/( W - Zf - t 2 (6) 

where w — wq + ie. (In the original inverse scattering method, before it was 
applied to general relativity, the poles were constants.) It is of interest to note 
that the pole trajectory is the result of setting the coefficients of any second 
order poles to zero, so the method is by construction only treating simple poles. 
For the Kasner metric, we note finally 

Now that the complete metric is given, it is of interest to count the free param- 
eters: the complex pole w (eq. the two complex parameters mo a (eq. ||) and 
the real scaling factor C2 (eq. pj), making for seven real parameters. 

We will also need the expression for the one real pole solution || (recall the 
superscript "(1)" indicates one pole): 

"> - f(& + %£u*) CD 



■ib 



r(l) = Cl f(0) 



with the auxiliary quantities 



m a = (moi|2w A| i m 02 \2w il 
L a — ymig n , m 2 g 2 2 ) 



~ IS— 1\ 



- 2 (0) , - 2 (0) 



(J, = Wq 



-Z±y/(wo-z) 2 -t 2 



We will see that these are real parts of the corresponding variables for g^ above, 
e.g. L a — > L a as e — > 0. Again we count the parameters and find only m 0a , C\ 
and Wq, four real parameters. 



3 The real limit 

3.1 General conclusions 

The procedure for taking the limit is not unique. For instance, let us take the 
expression for the 2x2 scattering matrix which generates the transformation, 
given by Belinskii and Zakharov fa] as 

x(A,M) = l + -^ + -^ (8) 

A — 11 A — LI 
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where Ri and R2 are 2x2 matrices and A is the spectral parameter. Now, if 
[X —>■ p, as e — > 0, we would obtain 

x(A,2,t) = l + . 



X — fi 



indistinguishable from the expression for one real pole. This seems to verify 
that the limiting procedure will, indeed, produce a simple real pole. However, 
one should note that the original expression (^) was intended for distinct poles 

Instead of taking the limit in m) we consider the behavior of the actual 
metric tensor as the complex conjugate poles wq ± ie approach the real pole wq 
by sending e — » (see fig. 0). 



Figure 2: Taking the limit e — > in w — wq ± ie. 



It is useful to discuss regions I (inside the light cone (wq — z) 2 < t 2 ) and 
II (outside the light cone (wq ~ z) 2 > t 2 ) separately. Refer to fig. [I] for the 
following discussion. 

i) Region I 

From equation (|J), region I will never have a real pole trajectory, but 
instead g approaches the background. Indeed, from that equation one 
finds 

|/i| 2 — > t 2 when e — > (region I) , 

which implies that g approaches the background. This in itself would 
suggest that we match the solution to the background in region I. On the 
other hand, 

gtt = ~9zz = f(z, t; e) 00 as 1/e 2 (region I) 

but as long as e is nonzero, one can still rescale / using the scaling constant 
C2 in equation (|§|) to allow for matching. 

ii) Region II 

We now expand the pole trajectory \x (eq. |^) in region II and find 

(j,(z, t) — ► jl(z, t) + in(z, t) e (region II) , 

where n{z, t) does not approach zero as e — * 0. Thus the complex conjugate 
pole trajectory does approach the corresponding trajectory for the real 
pole solution. 
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Now for the metric itself in region II. Special attention must be paid to 
the parameters; as was noted earlier, simple real pole solutions have only 
four real parameters while the complex pole solution has seven. We know 
the pole w becomes real. The two complex parameters mo a must become 
real in the limit, which accounts for the difference. An obvious ansatz for 
m Qa is 

m Qa = m 0a + ie p 

i.e. to let the complex parameters approach the corresponding real pa- 
rameters as some power p of e. With this ansatz, it can be seen from 
equation (||) that 

V -> as ae 2 + be 2p (region II) 

where a and b do not approach zero. 

We find two different cases: for p < 1, the second term in (g) vanishes, 
since T> approaches zero slowly - we are left with only the background in 
region II (as well as region I). On the other hand, for p > 1 we do obtain 
a new solution in region II. 

These results differ somewhat from earlier assertions 0, which state that 
the parameters must be proportional to each other, otherwise a finite pertur- 
bation cannot be found in the limit. This is again related to the issue of the 
number of free parameters, which is reduced when the parameters are taken to 
be proportional. Complex conjugate parameters can never be proportional, so 
if one takes this claim literally the limiting procedure can only be applied to 
the case of two real poles. However, in this argument the behavior of the pa- 
rameters niQa was not considered, and our scheme relies on the fact that these 
parameters approach reality with the poles, so it is free of this restriction. In 
[ pi] it is shown that a double pole can be obtained from two coalescing poles if 
once again the parameters are proportional, but the paper does not study the 
question whether single poles can be formed or not. Neither [Q nor jl6| treated 
the possibility of taking the real limit of complex constants. 

Now we turn to discuss the limit in detail. The general non-diagonal case 
becomes rather cumbersome, so we perform the limit numerically. As an in- 
structive example, we begin by studying the diagonal case analytically. 



3.2 The diagonal case 

A diagonal perturbation is obtained by taking, for example, mQ2 = 0. For a 
pair of distinct complex conjugate poles, we find 

t 2 (o) MM (o) 

511 = — g n , #22 - 922 ! 

f = C2 Im^u) 2 f (0) ( 9 ) 

\2w^i\ 4s (t 2 — \[i\ 2 ) 2 (t 2 — /i 2 )(i 2 — M 2 ) 
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which is known as the Einstein-Rosen metricHj. At this point, one notes that 
the exact same expression applies to the case of two real poles if /ii and /j,2 are 
substituted for /i and /2, respectively. If we take the diagonal perturbation for 
one real pole, on the other hand, we find from equation (Q): 

- 1 n - M J 

3n — Ul^ 11 ' 922 -— 922 ) 



2 (0) 

\2w fi\ 2s (t 2 - M 2 )Vi ' 



By inspection, we cannot obtain this metric from (^) by e — > 0. On the other 
hand, the two real pole metric has exactly the form M). Therefore, if the limit 
is taken at this stage, one does not obtain a simple real pole but two confluent 
real poles. This conclusion agrees with an older study Jl7| , which was limited to 
the present special case of a diagonal perturbation. The solution is, in general, 
not the two real pole solution itself since there is again ambiguity in the way the 
parameters coalesce. In addition, the expression for two real poles was intended 
for distinct poles, but this solution clearly exists as a limit. This is what we 
mean by "a pole on a pole" . 



3.3 The non-diagonal case 

We perform the limit numerically in the non-diagonal case. Details are avail- 
able elsewhere [[l8] , but the salient points are summarized here. We find that 
conclusions from the diagonal case carry through to the non-diagonal case. The 
metric for complex conjugate poles approaches the background everywhere for 
p < 1, but approaches a finite perturbation for p > 1, just as asserted on ana- 
lytical grounds earlier by studying T>. An example of this is given in fig. |3| for 
p = 2. The two complex pole solution clearly does not approach the one real 
pole solution. 

It appears that the two complex pole solution approaches the two real pole 
solution in the limit instead. However, the particular parameter dependence 
(given in the caption) has been chosen such that this would happen; in general 
the limit is different, although it is true that the limit procedure shown in fig. 
H may be thought of as the most natural limit for this comparison. All this is 
in agreement with the conclusion in the diagonal case. 

It is easy to see numerically that the metric coefficient / does approach 
infinity as e~ 2 . Belinskii and Francaviglia point out that if the scaling con- 
stant C2 is replaced with (essentially) e 2 , we can match / to the background 
on the light cone @. This would seem to allow for a finite / in the limit, since 
lim e _,o e2 / e2 = 1- Unfortunately, when e — 0, this would also mean = 0. 
Since there is only one scaling constant, the metric would then become degen- 
erate in region II, where / does not approach infinity. 

We conclude that the limit of two complex conjugate poles does produce a 
solution, but it is not the one real pole solution. Furthermore, this "pole on a 
pole" solution also suffers from the deficiency of being infinite inside the light 
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cone. Thus, the problem of finding a coordinate transformation which removes 
the infinity in / remains for the one real pole solution. 

4 Nonexistence of improved transformations 

Instead of the limit procedure above, one could attempt to improve the coor- 
dinate transformation for a simple real pole which has been used by several 
authors ||, ^] . A nonsingular transformation, which still does the job of render- 
ing / finite, would eliminate the null fluid. 

We here briefly summarize the results from Curir and one of the authors B. 
In terms of null coordinates 

u = t + z — Wq , v = t — z + Wq , 

the metric can be written as 

ds 2 = f dudv ~ g ab dx a dx b , a, 6=1,2. 
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When approaching the light cone u = from the perturbed region II, / goes 
to infinity as \j\f—u (v = can be treated in a similar way). By changing 
coordinate from u to 

, ciQo / 



V2 V 

where Qo is the value of Q on the light cone, the metric changes to 

ds 2 = /' du'dv - g ab dx a dx b , 

where /' is finite on the light cone and matches with /(°) from the unperturbed 
region I: f(u' = 0) = /<°)(u = 0). (By defining v! = u in the unperturbed 
region I, one sees that the light cone still is given by u' = and that v! is 
continuous across the light cone.) All metric components are now continuous, 
but their first derivatives suffer step function jumps across the light cone, giving 
rise to delta- functions both in the Weyl tensor and the Ricci tensor. 

So far points with the same value of v! and v have been identified along the 
light cone. Could one make the discontinuities disappear by identifying them in 
some other way? We have the following result: 

Proposition 1 There is no (topology-preserving) coordinate transformation of 
the one real pole solution (eq. |?|) which both makes f everywhere finite and 
removes the 8 -functions in the Ricci tensor. 

We proceed to show this. Clearly the light cone is given by u' = 0, so the 
remaining freedom is to identify different v with each other. But since the 
2-dimensional matrix g a t is continuous across the light cone one finds 



g a bi{u' = 0, vj) = g a bii{u' = 0, v n ) , 



implying that vj — vn. New coordinates u and v (we choose to perform a 
transformation in region II only, since this is general enough) must then satisfy 

u(0, v) = u' = , v(0, v) = v , 

on the light cone. Allowed transformations can then be written as 

u = u + F(u, v) , v = v + G(u, v) , 

where F and G are zero on the light cone but otherwise arbitrary functions. The 
requirement that the metric remains continuous then gives that F and G must 
approach zero faster than u. On the other hand, to avoid the discontinuities 
giving rise to the delta-function in the Ricci-tensor, G must approach zero as u 
or F must go as yj—u. Hence the we cannot avoid a distributional- valued Ricci 
tensor by this transformation. 

We do not consider transformations which induce topology change, such as 
rotating the light cone before identifying points. Such transformations would 
create completely different spacetimes than the ones we consider here. 

Our conclusion is that the null fluid is a property of the one real pole solution 
which cannot be removed. 
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5 Conclusion 



The original idea of eliminating the null fluid in the one real pole solution 
through a limit of two complex conjugate poles w = wo±ie is actually impossible 
to realize. We reach a different solution, "a pole on a pole" instead of a simple 
real pole solution. Also, we found that care must be exercised in the taking of 
the limit since the parameters must turn real. This is evidenced by the fact that 
in one case (what we call p < 1) the inverse scattering transformation reduces 
to the identity transformation. 

Additionally, we show that there is no transformation which accomplishes 
the goal of removing the metric infinities while keeping the null fluid away. 
Therefore, barring new topologies, the simple real pole solution really does con- 
tain a null fluid with negative energy density. We find this property unphysical 
in a completely classical solution, one would therefore reject these solutions. 
However, a more liberal interpretation would be that negative energy on the 
wave fronts is a characteristic of cosmological waves, and that this is a sign that 
quantum effects in the early universe are required to produce cosmological shock 
waves if such waves are produced at all. 

Since the metric in region I approaches the background under the limiting 
procedure, we found it natural to match region II to the background. It makes 
the wave interpretation more meaningful since the metric will be asymptotically 
Kasner on both sides of the wave front. It is also common to do so|l[. However, 
it is of course possible to match the solution in region II to some metric other 
than the background. See e.g. j^, for treatments where solutions generated 
by the inverse scattering technique from Bianchi models are matched to metrics 
other than the background. 

Higher-order (i.e. non-simple) poles in the scattering matrix could also gen- 
erate viable solutions, but that would be a different story altogether. 
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